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We apply time-dependent density-functional theory to study many-electron transport in 
Aharonov-Bohm interferometers in a non-equilibrium situation. The conductance properties in 
the system are complex and depend on the enclosed magnetic flux in the interferometer, the number 
of interacting particles, and the mutual distance of the transport channels at the points of encounter. 
Generally, the electron-electron interactions do not suppress the visibility of Aharonov-Bohm os- 
cillations if the interchannel distance - determined by the positioning of the incompressible strips 
through the external magnetic field - is optimized. However, the interactions also impose an inter- 
esting Aharonov-Bohm phase shift with channel distances below or above the optimal one. This 
effect is combined with suppressed oscillation amplitudes. We analyze these effects within different 
approximations for the exchange-correlation potential in time-dependent density-functional theory. 

PACS numbers: 73.63.Kv, 07.60.Ly, 31.15.Ge, 73.43. Cd 



I. INTRODUCTION 

Integer and fractional quantum Hall efFectsi together 
with the Aharonov-Bohm (AB) effecti^ enable the design 
of interferomctry devices that can be controlled by an ex- 
ternal magnetic field alone. Significant experimental ef- 
forts in two-dimensional (2D) systems have been made^"— 
to exploit the AB effect in the control of electron dynam- 
ics in quantum Hall interferometers. The complexity of 
these systems arises from the fact that not only the quan- 
tum mechanical phase of the electrons, but their path 
itself, is affected by the magnetic flux. Hence, the elec- 
tron transport in AB interferometers also provides a the- 
oretical challenge that should be approached beyond the 
single-particle edge-state scheme^ that does not take the 
path-dependence completely into account 

The characteristics of the measured AB oscillations 
have been analyzed with edge-channel experiments^iii 
and simulations In a recent work of some of the present 
authorsi^ the picture was made more complete by com- 
bining full electrostaticsi^ of a AB interferometer de- 
vice^ with a time-propagation of noninteracting excess 
electrons in the obtained incompressible strips in non- 
equilibrium. Clear AB oscillations were found as a func- 
tion of the magnetic flux, and the visibility (amplitude) 
of the oscillations was found to largely depend on the 
width and position of the current channels. 

In this work we extend the earlier single-electron 
studyi^ to the many-electron regime in the frame- 
work of time-dependent density-functional theor y ^^1^^ 
(TDDFT). We consider different (albeit adiabatic) ap- 
proximations for the exchange-correlation functional in- 
cluding the exact exchange (EXX), Hartree, and local- 
density approximation (LDA). We find that in the high- 
visibility regime, corresponding to an optimal positioning 
of the incompressible strips (and thus the current chan- 
nels), the electron-electron interactions do not suppress 



the AB oscillations. This is in accordance with a previ- 
ous TDDFT study on semiconductor quantum rings 
However, the interactions introduce an AB phase shift 
whose magnitude is an increasing function of the elec- 
tron number varied between one and ten. In addition, 
the interactions suppress the AB amplitudes if the cur- 
rent channels arc either too close or too far apart from 
each other. The induced phase shift could be measured 
in an experiment, which would eventually quantify the 
actual role of electronic interactions in a real setup. 



II. MODEL AND METHODS 

The transport simulation is carried out by injecting 
many-electron wave packets into a static 2D potential. 
Figure [T] shows a schematic picture of the system. The 
shape of the model potential is based on the experimen- 
tal GaAs/AlGaAs device by Camino et alA and electro- 
static calculations performed in Ref. [Tsl . In particular, 
the potential corresponds to incompressible strips with 
an integer filling factor (here v — 2) that carry the non- 
equilibrium electric current. The model potential along- 
side the thick curves in Fig. [T] has a Gaussian cross sec- 



tion, Vci 



-Vbe 1'^ , where s is the coordinate per- 



pendicular to the channel, Vb = 80 (Hartree atomic units 
used throughout) is the channel depth set to prevent the 
electrons from escaping the tracks, and c = 0.3 is the 
track width parameter. A part of the potential is visual- 
ized as a mesh in the upper part of Fig. [TJ A Gaussian 
cross section is a convenient choice for two reasons. First, 
it is a good approximation for the magnetic (parabolic) 
confinement at the bottom of the channel. Secondly, the 
upper part of the potential allows "leaking" of the elec- 
tron flow in consistency with the experimental situation. 
It is important to note that in the time-propagation we 
deal with the excess electron density on top of the Fermi 
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FIG. 1: (Color online) Schematic figure of the transport setup. 
The upper mesh represents the cross section of the model 
channel potential along the blue (black) thick lines. The 
lower mesh shows the initial wave packet used in the time- 
propagation. 

background that is assumed to be static. 

We also apply a linear ramp potential to accelerate 
the initial state, thus mimicking the source-drain voltage. 
The initial state in the calculation is an iV-electron wave 
packet at the lower-left corner of the system (see the lower 
mesh in Fig. [T]). It is calculated by solving the ground- 
state problem for a part of the input channel (confined 
by step functions) that contains N Coulomb-interacting 
particlesji^ 

Our transport scheme described above leads to in- 
evitable non-chiral transport, i.e., the direction of the 
current is the same on the both sides of the interferome- 
ter. This is in accordance with non- equilibrium transport 
experiments)^ where an external current is imposed and 
a Hall (electrochemical) potential difference develops at 
two opposite edges with the same sloped Such poten- 
tial drops in the incompressible strip have been found in 
several experiments)^"— In the conventional equilibrium 
case, however, the current would be chiral as confirmed 
in several studies. 

In the time propagation of electrons we apply 
TDDFTi^iii on a 2D real-space grid. We use the fourth- 
order Taylor expansion of the time evolution opera- 
tor to the Kohn-Sham states and perform the calcula- 
tion with the OCTOPUS code packagei^ The exchange- 
correlation potential in TDDFT is treated within adia- 
batic LDA, and in the two-particle case also with the 
(adiabatic) exact exchange approximation. In compar- 
ison, we have applied the bare Hartree approximation 
(including only classical electrostatic interactions and 
neglecting exchange and correlation) and independent- 
electron approximation (neglecting interactions com- 
pletely). 

We consider the conductance of the system as a func- 
tion of the uniform and perpendicular magnetic-field 
strength in the units of the enclosed flux quanta, <&/$0j 
where $ = i? x area and $o = h/e. Another changing 
variable is the distance between the points of encounter 
(interchanncl distance) d (see Fig. [T|) and the number of 



injected particles TV = 2, 6, or 10. We consider only 
spin-compensated (zero total spin) systems. We point 
out that in an experimental situation d would be essen- 
tially controlled either with the gates, or with the mag- 
netic field that affects the position of the incompressible 
strips. However, in our transport model we allow the 
magnetic field affect only the flux but not the distribu- 
tion of the channels; this is justified by the fact that with 
realistic parameter values our changes in the flux would 
correspond to only minor changes in the magnetic field 
with respect to the effect needed to change the spatial 
distribution of the incompressible strips. 

The relative conductance is estimated as the number 
of electrons found in the top-right corner of the system 
as a function of timcj^ Since the simulation box is fi- 
nite, wc monitor the electron density only until back- 
scattering effects from the borders of the box appear. 
The measured conductance parameter is the normalized 
probability density found in the upper right corner of the 
system. A more detailed discussion of the model system 
can be found in Ref . [l^ 



III. RESULTS 

First wc consider electron transport in a two-electron 
system. The initial state is injected in the interferometer 
using the linear ramp potential, and the relative conduc- 
tance is assessed as described in the previous section. The 
panels in Fig. [5] show the conductance for different inter- 
channel distances, respectively, and different approxima- 
tions for the electron-electron interactions are compared. 

In all cases we find distinctive AB oscillations with 
the correct flux quantization in the noninteracting limit. 
Hence, the direction of the electron current (to the left 
or right drain) is completely controlled by the flux. The 
maximum amplitude (visibility) of the AB oscillations 
strongly depends on the interchanncl distance. Regard- 
less of the interactions, the maximum visibility is found 
at c? ~ 0.4. A detailed look in the electron densities show 
that in this case the partitioning of the wave packet in 
the first point of encounter (see Fig. [Ij is close to even, 
i.e., one electron remains on the left track, but the other 
one is transported to the right track. Consequently, the 
interference at the second point of encounter is as com- 
plete as possible. 

For two electrons the interactions seem to have a minor 
effect on transport characteristics. They only induce a 
small phase shift in the AB oscillations. The direction 
of the phase shift depends on d: at d < 0.4 (below the 
maximum visibility) the oscillations are shifted to higher 
<I>/<I>o due to interactions, and at c? > 0.4 the situation 
is opposite. The phase shifts are analyzed in more detail 
below within Fig. [S] 

Moreover, as seen in Fig. [2l the EXX and LDA re- 
sults are generally very close to each other. As the elec- 
tronic correlations are expected to be relatively small, 
the obtained agreement validates the accuracy of LDA 
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FIG. 2: (Color online) Estimated conductance in two-electron 
Aharonov-Bohm interferometers with different interchannel 
distances and different approximations within time-dependent 
density-functional theory (independent electrons, adiabatic 
local-density approximation, adiabatic exact-exchange calcu- 
lation). The maximum visibility is obtained with d ~ 0.4. 



(as in EXX the exchange part is exact). In this respect, 
our dynamical results are in accordance with previous 
density-functional studies where the (surprisingly) high 
accuracy of LDA results have been pointed out - even 
in complicated situations such as quantum Hall droplets 
with fractional filling factors such as = 5/2 (Refs.ii- 
13) and z/ = 1/3 (Ref.Hi). 

In Fig. [3] wc plot the same data as in Fig. [5] but sepa- 
rate the results in terms of the interaction scheme. The 
homogeneity of the AB oscillations in the noninteracting 
case (a) is very clear, so that only the amplitude changes 
as a function of d. As discussed above, LDA and EXX re- 
sults are similar and both show considerable phase shifts 
as a function of d. It should be noted that the physical 
mechanism behind the phase shift is not in the different 
partitioning in the presence of interactions as this would 
mainly affect the amplitudes. Instead, the interactions 
enhance the spreading of the wave packet due to Coulomb 
repulsion in an asymmetric fashion, i.e., the spreading is 
relatively more prominent on the higher-density side of 
the interferometer. As the spreading is also radial, there 
is an asymmetry in the capturing of the phase that leads 
to a shift in the AB oscillation phase. 

In Fig. m we plot the conductance as a function of the 
flux for a ten-electron system. The results with distinc- 
tive AB oscillations are qualitatively similar to the two- 
electron situation considered above. However, there are 
some important differences. First, the maximum am- 
plitude is obtained with a smaller interchannel distance 
{d ^ 0.4) than in the two-electron case. This is due to 




FIG. 3: (Color online) Same two-electron data as in Fig.[2]but 
plotted separately for different approximations within time- 
dependent density-functional theory. 



the maximum partitioning (50% vs. 50%) to the left and 
right channels with this d, so that the underlying mech- 
anism is similar to the two-electron system. Secondly, 
we find that apart from the d = 0.4 case the ampli- 
tude is significantly dampened by the electron-electron 
interactions. Thirdly, as shown in Fig. [SJ also the AB 
phase is affected by the interactions more than in the 
two-electron system. This can be understood by the fact 
that the density cloud is larger, so that its parts capture 
different phases. Consequently, the overall amplitude is 
partly smeared out. It is noteworthy that in the case of 
d = 0.6 the AB oscillations are almost completely damp- 
ened in the LDA solution [see Fig.lSJc)]. 

The ten-electron system in Figs. U and [5] does not 
enable numerically stable EXX solutions. However, we 
compare the adiabatic LDA calculations with the Hartree 
results, where the exchange-correlation effects are omit- 
ted. The results in the high- amplitude regime are simi- 
lar, and also generally the deviations are relatively mod- 
est. Overall, the high-amplitude (high- visibility) regime 
is not considerably affected by the electron-electron in- 
teractions. Therefore, we might expect that in a real ex- 
perimental setup the AB oscillations arc always present 
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FIG. 4: (Color online) Estimated conductance in ten-electron 
Aharonov-Bohm interferometers with different interchannel 
distances and different approximations within time-dependent 
density-functional theory (independent electrons, adiabatic 
local-density approximation, Hartree approximation). 



if the partitioning is optimal, which is eventually deter- 
mined by the location of the incompressible strips con- 
trolled by the initial magnetic field. Moreover, the AB 
phase can be controlled by changing the positioning of 
the incompressible strips; however, the phase change is 
achieved only with the price of suppressing the visibility. 

In Fig. iniwe present a summarized visualization of the 
AB oscillation visibility with different number of elec- 
trons propagating in the interferometer, (a) N = 1, (b) 
iV = 6, and (c) N = 10. For the six- and ten-electron re- 
sults we use the adiabatic LDA. Our results demonstrate 
that, as discussed above, the overall visibility is not re- 
duced as a function of N. On the other hand, the phase 
shift around the "optimal" interchannel distance - that 
is always located around d = 0.4 . . . 0.5 - is proportional 
to N. Moreover, unless the interchannel distance is close 
to the optimal value, or below that, the interactions re- 
duce the AB oscillation amplitudes. This is visualized 
in Fig. [7] where we plot the maximum amplitude in the 
whole range <I>/<I>o — . . . 3. We find that in the inter- 
acting = 10 case the amplitude decays close to zero 
at d = 0.6, whereas at that distance the single-electron 
result - and the A^ = 10 without electron-electron inter- 
actions - still have a 20% visibility. 



IV. SUMMARY AND DISCUSSION 

In summary, we have applied a single-electron propaga- 
tion scheme as well as time-dependent density functional 
theory to study electron transport in realistic Aharanov- 
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FIG. 5: (Color online) Same ten-electron data as in Fig.|4]but 
now plotted separately for different approximations within 
time-dependent density-functional theory. 



Bohm interferometers in the integer filling factor regime. 
The device geometry has been obtained from experimen- 
tal boundary conditions that have been applied in a self- 
consistent Poisson solver within the Thomas-Fermi ap- 
proximation in order to obtain the positions of the in- 
compressible strips. Here we have assumed, in accor- 
dance with previous experiments and calculations, that 
the electric current is carried along the incompressible 
strips in a non-chiral fashion (as the system is in non- 
equilibrium). The current channels have been modeled 
with a Gaussian tubes where electronic wave packets have 
been injected. 

The time-evolution of excess electron density leads to 
distinctive Aharonov-Bohm oscillations as a function of 
magnetic flux. The nature of the oscillations depend on 
the electron-electron interactions handled within (adia- 
batic) exact-exchange, local-density approximation, and 
Hartree approximation, respectively, in time-dependent 
density-functional theory. We have found that the local- 
density approximation is accurate in comparison with the 
exact-exchange approach. 

The electron-electron interactions do not affect the 
maximum visibility even if the number of propagat- 
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FIG. 6: (Color online) Aharonov-Bohm oscillations (color 
scale) for one-, six-, and ten-electron systems, respectively. 
Clear phase shifts can be in many-electron systems treated 
here within adiabatic local-density approximation. 



ing electrons is increased to ten. This is in agreement 
with a previous time-dependent density-functional study 
on many-electron semiconductor quantum rings How- 
ever, we have found that around the optimal interchan- 
nel distance the interactions induce a phase shift in the 
Aharonov-Bohm oscillations, and the strength of the shift 
is proportional to the number of electrons. If the inter- 
channcl distance is far from optimal the interactions may 
dampen the Aharonov-Bohm oscillations completely. 

Our results confirm that in a corresponding experi- 
mental setup the Aharonov-Bohm oscillations should be 
clearly visible, especially if the incompressible strips are 
sufficiently narrow. Furthermore, possible phase shifts 
in the oscillations in different magnetic- field ranges, i.e., 
with different positions and widths of the incompressible 
strips, could reveal the importance of interaction effects 
in a real setup. As a signature of such a dependency 
on the interchannel distance we would like to mention 
the early experiments performed in Mach-Zehnder ge- 
ometryfii There the visibility of Aharonov-Bohm oscilla- 
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FIG. 7: (Color online) Maximum amplitude of the Aharonov- 
Bohm oscillations (from minimum to maximum conductance) 
in the range "l?/<&o = ... 3 in the case of A*' = 1 (solid line), 
A" = 10 caluculated with the adiabatic local-density approxi- 
mation (dashed line) and without interactions (circles). 

tions is suppressed depending on the current excitation, 
and this phenomena is attributed to the energy depen- 
dence of the transmission probabilities. We hope that the 
present study encourages further experimental efforts in 
this direction. 
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